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We analyze thermally induced spin and charge transport in HgTe/CdTe quantum wells on the 
basis of the numerical non-equilibrium Green's function technique in the linear response regime. In 
the topologically non-trivial regime, we find a clear signature of the gap of the edge states due to 
their finite overlap from opposite sample boundaries - both in the charge Seebeck and spin Nernst 
signal. We are able to fully understand the physical origin of the thermoelectric transport signatures 
of edge and bulk states based on simple analytical models. Interestingly, we derive that the spin 
Nernst signal is related to the spin Hall conductance by a Mott-like relation which is exact to all 
orders in the temperature difference between the warm and the cold reservoir. 



I. INTRODUCTION 

Thermoelectric transport coefficients define the effi- 
ciency of a system to generate an electrical power from 
a temperature gradient pQ. The most established ther- 
moelectric phenomenon is the Seebeck effect [2J [3J, in 
which a current (closed boundary conditions), or a bias 
(open boundary conditions) is induced from a tempera- 
ture difference held between two reservoirs of a junction. 
The transverse Seebeck coefficient, or Nernst coefficient, 
refers to the alternative situation where the thermally in- 
duced current (bias) flows in the direction transverse to 
both the temperature gradient and the applied magnetic 
field g]. 

Thermoelectric effects have major consequences in 
terms of technological impact and scientific understand- 
ing. On the one hand, these effects offer interesting appli- 
cations based on heat-voltage conversion: thermometry, 
refrigeration, power generation [6]. On the other hand, 
thermoelectric coefficients combine information from en- 
ergy and charge flows at quasi-equilibrium. Furthermore, 
they are more sensitive to the details of the density of 
states than electrical conductance [7H9]- Both aspects 
make them a powerful tool to probe the system dynam- 
ics [ID]. 

During the last two decades, there have been consider- 
able technological advances in low-temperature nanoscale 
physics. This allows precise measurements of thermo- 
electric transport signals, obtained in various systems 
like bismuth [TT] , superconductors [TU H3J , carbon-based 
structures [Uj, or molecular junctions [T5] . 

The recent alliance of spintronics and thermoelectric 
transport brings up a spin analog of Seebeck and Nernst 
effects (see Ref. [H] for a short review). Especially in sys- 
tems with strong spin-orbit interactions, a temperature 
gradient can generate a transverse spin current (or bias) 
even in the absence of an applied magnetic field. This 
can lead to the anomalous Nernst effect (in the case of 
ferromagnetic systems) |17ff20] or the spin Nernst effect 
(in the time-reversal symmetric situation) [2 H - E3] . 

Systems with strong spin-orbit interactions have been 
extensively studied in condensed matter physics espe- 
cially since the prediction of the spin Hall effect [24Tf2"T] 



which allows for an all-electrical manipulation of spin. 
The spin Hall effect generates a transverse spin accumu- 
lation as a response to a longitudinal applied electric field. 
Spin-orbit interactions have several origins which distin- 
guish the different types of phenomena, for instance, an 
extrinsic spin Hall effect can emerge from the spin-orbit 
dependent scattering on impurities or defects [Ml 12H1 
I31j . On the other hand, bulk or structure inversion 
asymmetries [32]-|34j give rise to an intrinsic spin Hall 
effect [2H1 (HI 135], which may be described in terms of 
an anomalous velocity or a spin-dependent classical force 
[36H38]. 

Recent experiments have demonstrated the existence of 
an intrinsic spin Hall effect in HgTe/CdTe quantum wells 
(QWs) [3J5] by the use of the quantum spin Hall effect as 
the detector. This novel electronic phase is characterized 
by an insulating bulk and protected metallic edge states. 
The emergence of the quantum spin Hall effect is due 
to strong spin-orbit coupling and other relativistic cor- 
rections, which reverse bands of opposite parities. The 
electrons obey a massive Dirac equation and the sign of 
the mass term enables us to distinguish the topological 
phases. The edge states consist of Kramers pairs moving 
in opposite direction at each boundary H0HI2] and time- 
reversal symmetry protects them from non-magnetic and 
elastic backscattering [53]. Thereby, these edge channels 
carry "dissipation-less" spin currents whose existence in 
HgTe QWs has been confirmed experimentally by mea- 
surements in multi-terminal devices [44] , 

Recently, topological insulators have been proposed 
as good materials for thermoelectric conversion [4"5lJ48| . 
The basic idea relies on the topological protection of ID 
edge states which prevents reduction of electrical trans- 
port in disordered systems. The authors consider narrow 
ribbons of quantum spin Hall insulators or 3D topological 
insulator with line dislocations. The aim is the enhance- 
ment of the contribution of edge states to the thermoelec- 
tric transport compared to the bulk modes. Hence, the 
analysis is restricted to a small energy range excluding 
the valence bands. Inelastic processes are taken into ac- 
count through the Boltzmann transport theory. The lat- 
ter is used to calculate the thermoelectric coefficients in 
a two-band model. Thus, the efficiency of these systems 
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to convert heat into electricity is based on the dominance 
of the edge modes on transport. 

In this work, we investigate the spin-dependent ther- 
moelectric transport in quantum spin Hall insulators 
based on HgTe/CdTe QWs in absence of magnetic fields. 
The behavior of the Seebeck coefficient and the spin 
Nernst signal is analyzed in a four-terminal cross-bar 
setup, as shown in Fig. [T] A thermal gradient between 
lateral leads induces a longitudinal electrical bias and a 
transverse spin current. Each of them can be used as a 
probe of the topological regime as well as finite size ef- 
fects of the quantum spin Hall insulator. We show that 
the oscillatory character of the Seebeck and spin Nernst 
coefficients in the bulk gap highlights the presence of the 
mini-gap - due to the finite overlap of the edge states 
from opposite sample boundaries. Furthermore, we de- 
scribe a qualitative relation between the type of particles 
in a given band and the magnitude of the spin Nernst sig- 
nal. This allows us to provide a natural explanation of 
the observed phenomena based on anomalous velocities 
and spin-dependent scattering off sample boundaries. 

The article is organized as follows. In Sec. II, we intro- 
duce the model Hamiltonian of the HgTe/CdTe QW and 
describe the formalism necessary to calculate the Seebeck 
and spin Nernst coefficients. The thermoelectric trans- 
port by the edge states - with a particular emphasis on 
finite size effects - is analyzed in Sec. Ill through the 
behavior of Seebeck and spin Nernst signals. In Sec. IV, 
we focus on the spin-dependent thermoelectric effect in- 
duced by the bulk states. We conclude in Sec. V and 
present details of the calculation in the appendices. 



II. MODEL 

In this section, we present the model Hamiltonian of 
the HgTe/CdTe QW and give the general expressions of 
the Seebeck and spin Nernst coefficients. 



A. Hamiltonian 

We consider a four-terminal cross-bar setup based on a 
HgTe/CdTe QW whose low-energy dynamics is described 
by the Bernevig-Hughes-Zhang (BHZ) 4-band model [42j 
33] . The Hamiltonian is written in the basis of the lowest 
QW subbands \E+), \H+), \E-), and \H-). Here, ± 
stands for two Kramers partners but in the following, we 
will refer to them as spin components, denoted by 
for brevity. The spin z-direction corresponds to the QW 
growth direction, which is [001]. The Hamiltonian can 
be written as 
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FIG. 1. (Color online) Four-terminal cross-bar setup based 
on HgTe/CdTe QWs used for thermoelectric transport. A lon- 
gitudinal temperature gradient AT is applied between reser- 
voirs a and b and generates a transverse spin current 1^ de- 
tected, for instance, in reservoir c. 



■ky, k± — k x ±ik y , and M(k) = M — Bk 2 



with k 2 = k 2 

The sign of the gap parameter M determines whether 
we are in the regime of a trivial insulator (M > 0) or 
a topological insulator (M < 0). Experimentally, M is 
tuned by changing the QW width. 

The term V m (r)r 2 describes an in-plane confinement 
potential, where t z is a Pauli matrix acting on the E/H 
space. By this kind of confinement we may ensure that 
outside of the sample, i.e. in vacuum, the parameter 
regime is topologically trivial, so that edge states, if 
present, will be confined. Calling the inside of the sample 
G, the limit V m (r) — > oo Vr e dG can be used to make 
all components of ip vanish at the sample boundary, the 
envelope function if) being the solution of the Dirac equa- 
tion based on the Hamiltonian Q. 

We mention in passing that this model can be ex- 
tended by a term breaking the structural inversion asym- 
metry (SIA) with a z-dependent potential. The resulting 
Rashba-like interaction connects the Kramers blocks of 
the Hamiltonian affecting the particles with opposite 
spin 
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with the Rashba coupling parameters i?oi $0, an d T [49j . 
We have analyzed that such a term will only quantita- 
tively affect all our results presented below. Therefore, 
we will not further consider effects due to SIA in this 
article. 

Figure [T] shows the four-terminal cross-bar setup we 
analyze. The central sample is connected to four semi- 
infinite leads: the reservoirs a and b are maintained re- 
spectively at warmer and colder temperature than the 
rest of the system - creating a longitudinal temperature 
gradient - while the transverse terminals c and d are used 
to probe spin currents. 

To model the setup and treat the thermoelectric trans- 
port properties, we employ the tight binding approach. 
Therefore, we discretize the continuum model iQ on a 
lattice of spacing a by the substitutions k 2 — > 3^(2 — 
2 cos kia) and ki - sin(fcia), where i is the index of the 
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tially overlap [50] . One of the consequences of this is the 
opening of a mini-gap, as shown in the inset of Fig. [2^b) . 



B. Landauer Biittiker formalism and 
thermoelectric coefficients 

The particle current l pa in the lead p with spin a is 
obtained by the Landauer-Buttikcr formula [51j 
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FIG. 2. (Color online) Subband dispersion relation for the 
leads (400 nm wide), in (a) the normal regime (M = O.lto) 
and (b) the inverted regime (M = —O.lto)- The inset shows 
the mini-gap in the dispersion caused by the overlap of edge 
states. The coloring highlights the transition from electron- 
like character (\E±) in red) to heavy-hole character (\H±) in 
blue) with full color for a 20% excess of either contribution. 
The lattice constant is a = 6.6 nm so to = 44 meV. The 
mini-gap width is 6.5 x 10" 3 to(= 0.28 meV). 



lattice site. The confinement potential is implemented 
by the lattice truncation in accordance with the geom- 
etry of the sample. Rewriting the trigonometric func- 
tions in terms of translation operators, this leads to a 
tight binding Hamiltonian which only contains nearest- 
neighbor hopping terms between the lattice sites (see 
Appendix [A] for details). The energies of the model 
are expressed as functions of the conduction band hop- 
ping parameter *o = — D ^~J* where the parameter val- 
ues of the HgTe/CdTe QW are taken as in typical ex- 
periments, i.e. A — 0.375 eVnm, B = —1.120 cVnm 2 , 
and D = —0.730 eVnm 2 . In the low-energy regime, the 
lattice constant is set sufficiently small compared to the 
Fermi wave length. Hence, for a = 6.6 nm, the energy 
unit is <o — 44 meV. The parameter M is chosen as 
\M\ = O.lto = 4.4 meV. In Fig. ^a) and (b), we show 
the subband dispersion relation for a HgTe/CdTe QW 
of width 400 nm, both in the normal insulator and the 
topological insulator regimes. In the latter case, finite 
size effects emerge on the edge states since they substan- 



tribution function, ks the Boltzmann constant, T p the 
temperature, and fj, p the chemical potential. The trans- 
mission probability T paq (E) — T paq(T i(E) from lead 
p with spin a to lead q can be evaluated using the non- 
equilibrium Green's function formalism (NEGF) [52H54] 



T P a, g *'(E) = Tr[r pCT G i 



Tgo-'G 



(p,a)^(q,a / ), (5) 



where T pa (E) = i(S pcr — E^) refers to projectors on ve- 
locity operators of the propagating modes, and T, pa (E) 
stands for the spin-dependent self-energy. The latter 



is defined by T, pa (E) = T pa (E + i0 + - #i ea ds) 



per 5 

where the matrix T pa connects the lead p, spin a to 
the adjacent sites of the sample. G R (E) = (G A (E)) J< = 
(E - -ffsampic - J2 P .a ^W)" 1 is tne retarded Green's func- 
tion. Further, i?i ca ds and -ff S ampie represent, respectively, 
the lattice Hamiltonians of the decoupled leads and the 
sample. Once the transmission probabilities are evalu- 
ated, the charge current I p — e(I p ^ + Ip^) and the spin 
current I p — (h/2)(I p ^ — I p ±) can be obtained. 

We consider a longitudinal temperature gradient be- 
tween the leads a and b by setting T a = T + AT, 
T b = T - AT, and T c = T d = T. In the linear response 
regime, the Seebeck coefficient reports the longitudinal 
voltage bias A/z = ^""^ generated by the temperature 
gradient AT under the condition of vanishing charge cur- 
rents (open boundary conditions). Upon Taylor expan- 
sion of the Fermi functions in AT and Ap,, the Seebeck 
coefficient can be written as 



S K = 



-A/x 



eAT 



1 fdEf (l-f )Ts E (E)(E-ri 
eT J dE f (l - f )TsE(E) 



with f = (e 



= ( P (E-n)/k B T_ 



(6) 

1) _1 the Fermi distribution func- 



tion at equilibrium. In the above equation, we defined the 
Seebeck transmission function Tse(E) — T 3 b + (T a . c + 
T a ,d)/2, where the summation over spins is implied. 

Due to the presence of intrinsic spin-orbit interaction 
in the sample, the longitudinal thermal gradient AT also 
yields a transverse spin current —1$) in the case of 
closed boundary conditions. The spin Nernst coefficient 
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is then defined as the ratio 



2AT 
1 



8ttA: b T 2 



d#/o(l- fo)T SN (E)(E - fx). (7) 



Here, we introduced the spin Nernst transmission func- 
tion Tsn(E) = AT c .b — AT ca , with the short- hand nota- 
tion AT Cib (£) = T ctM + Tcf^i - T ciM - T ci:hi . 

Interestingly : the Mott relation provides information 
about the (spin-)thermotransport coefficients on the ba- 
sis of the energy dependence of the (spin-)conductance 
|55j . In the low temperature limit, one can derive 



7r 2 fc|T dlnG xx (E) 



3e 



dE 



(8) 



where G XX (E,T = 0) = \T S e(E) is the longitudinal 
conductance for zero temperature. Equation ^ is valid 
if /csT is large compared to the scale on which Tse{E) 
varies. A numerical analysis in 56j claims that Q can 
be valid even if Tse(E) varies more rapidly, as long as 

An analogous relation exists between the spin Nernst 
signal and the spin Hall conductance. From the Sommer- 
feld expansion of the transmission function in Eq. ([7]), 
one obtains the following Mott-like formula 



2n 2 k 2 B T dG sH (E) 



3< 



dE 



(9) 



B =M (T=0) 



with G sH (E,T = 0) = ^Tsn(E), the spin Hall conduc- 
tance at zero temperature. In Appendix [5J we demon- 
strate that this relation can be extended to finite temper- 
ature by defining a smoothed function Tsn(E) (see Eq. 
(B5)) that depends on the temperatures in the leads. As 
a result, we find an exact Mott-like formula for the spin 
Nernst coefficient 



7Tfc|T dt SN (E) 
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dE 



(10) 



Since Tsn(E) shows a highly oscillating behavior, the 
above equation simplifies the interpretation of the spin 
Nernst signal in terms of transmission f unc tions because 
of the smoothing of Tsn(E). Equation (10) is one of the 
key results of our paper. 



III. THERMOELECTRIC TRANSPORT 
CARRIED BY THE EDGE STATES 

In this section, we present the numerical results of the 
spin Nernst and Seebeck coefficients for an energy regime 
within the bulk gap. When the HgTe/CdTe QW is in a 
topologically trivial phase, there is no sub-gap transport 
through the system. The transmission functions Ts e and 
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FIG. 3. (Color online) Four-terminal cross-bar setup based 
on a HgTe/CdTe QW in the inverted regime. When an elec- 
trical bias is imposed between the longitudinal leads a and b, 
one edge channel carries electrons with spin up from lead a to 
lead c (red solid line) and one edge channel carries electrons 
with spin down from lead b to lead c (green solid line). This 
gives rise to a spin current in lead c and, at zero temperature, 
to a quantized spin Hall conductance. 



Ts n are then zero, and from Eqs. ([6| and Q, it follows 
that there are no thermoelectric signals. On the contrary, 
the HgTe/CdTe QW in a non trivial phase hosts edge 
states in the bulk insulating gap. These modes carry 
electrons with opposite spins in opposite directions. 

With respect to the geometry of the setup, spin and 
electrical currents are induced and flow respectively in 
transverse and longitudinal leads, as depicted in Fig. [3] 
However, the finite width of the system implies an overlap 
of the edge states meaning that backscattering processes 
can occur. 

We first investigate the behavior of the spin Nernst 
signal N s and the associated transmission function Tsn- 
The results are presented in Figs. Qa) and (b). While the 
chemical potential is in the bulk gap, the spin transport 
is mediated by the edge channels so that the transmission 
function is simply given by 

Tsn = (T ct>b - T cl . b ) - (T cT , - T 4>a ) = -2. (11) 

Evidently, as the chemical potential reaches the bound- 
ary of the mini-gap, the number of propagating states 
drops to zero and transport breaks down. This results in 
a peak of the transmission function Tsn- Consequently, 
the spin Nernst coefficient is zero in the bulk gap except 
when the chemical potential reaches the boundary of the 
mini-gap. Because of the Mott-like relation (10), a sym- 
metric function Tsn(E) must result in an antisymmet- 
ric function N S (E). Therefore, N s exhibits an approxi- 
mately antisymmetric peak centered at the maximum of 
the transmission peak. The confinement of the QW im- 
plies an energy shift in the band dispersion. Therefore, 
the boundaries of the bulk gap are not exactly at energy 
\M\ = O.lio, as we can see in Figs. [2] and [4] In Fig. 
|4j the gap and minigap positions of a finite system are 
indicated by vertical lines. Dotted vertical lines are used 
for the inverted regime and dashed lines for the normal 
regime. Interestingly, one observes that the merging of 
the edge state to the conduction band causes Tsn to van- 
ish already before the first bulk mode appears. Where fi 
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FIG. 4. (Color online) (a) Spin Nernst transmission function 
as a function of the chemical potential when the sample is in 
the normal (dashed blue line) or in the inverted (solid red line) 
regime. The dotted vertical lines indicate the bulk gap and 
minigap positions in a finite system for the inverted regime, 
while the dashed vertical lines indicate the gap in the normal 
regime, (b) Spin Nernst signal N a /kB in a system in the 
inverted regime at T = 2 K (black solid line), T = 6 K (green 
dashed line), and T = 8 K (orange dotted line). 

lies between the right dashed and dotted vertical lines, 
a finite Ts jv reappears due to the formation of the first 
bulk state, in the same subband as the edge state. Out- 
side the gap indicated by the vertical lines, bulk states 
start to participate to the spin transport resulting in ad- 
ditional oscillations in Tsn as a function of fi. They 
transform into peaks of the spin Nernst coefficient whose 
existence is understood with the same arguments as for 
the mini-gap peak. Especially at positive chemical po- 
tential, the magnitude of the peak is comparable to that 
of the mini-gap and allows to mark the position where 
the edge states merge. 

In Fig. ^b), we show the behavior of the spin Nernst 
coefficient for different temperatures. As fc^T increases, 
the position of the peaks is slightly shifted to lower en- 
ergy. The magnitude tends to decrease and the peak 
width is broadened. Up to T = 6 K, the spin Nernst 
signal goes to zero between the peak that specify the po- 
sition of the mini-gap and the edge state merging peak. 
Beyond this temperature, N s is smoothed out, so that it 
can not probe the edge state signal. 

We now turn to the analysis of the transmission func- 
tion Tse and the Seebeck coefficient S e as a function of 
energy. The results are presented in Fig. [5] Inside the 
bulk gap, the transmission function Tse is constant but 
goes to zero when the chemical potential is in the mini- 
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FIG. 5. (Color online) Seebeck coefficient 5 e [fcs/e] (red solid 
line) and scaled transmission function Tse/ 50 (black dashed 
line) as a function of the chemical potential at T = 2 K. The 
mini-gap appears as an antisymmetric peak. 

gap. This feature leads to an approximately antisymmet- 
ric peak in the behavior of S e , which provides information 
on the presence and the position of the mini-gap in the 
spectrum. The boundary of the bulk gap manifests itself 
as the step of the transmission function and transforms 
as a narrow peak in S e . 

The transmission function Tse exhibits a smoothed 
staircase behavior whose steps coincide with the opening 
(at positive energy) or the closing (at negative energy) 
of conducting channels. This behavior transforms into a 
series of peaks in S e . However, as the chemical poten- 
tial increases, the magnitude of the peaks reduces. The 
reason is that the considered setup possesses four termi- 
nals that all exhibit an increasing number of propagating 
modes with increasing /i. Thus, inter-mode scattering is 
more and more likely to happen. Then, the staircase be- 
havior of Tse diminishes and transforms into oscillations. 

We close this section with a remark on the average sign 
of the Seebeck coefficient S e ■ It is positive in the conduc- 
tion band and negative in the valence band which reflects 
the sign of the corresponding excitations in a given band. 

IV. SPIN NERNST EFFECT INDUCED BY 
BULK STATES 

A spatial dependence of model parameters, like, for in- 
stance, an in-plane electrostatic potential Vo or the mass 
confinement potential V m , can generate a transverse spin 
current resulting in a spin Hall signal for the metallic 
bulk states |49j . This phenomenon has been previously 
analyzed in Refs. [5 7-53] in the context of charge and 
spin transport properties at interfaces between metals 
and quantum spin Hall systems. As already mentioned 
above, the spin Hall conductance gives rise to the spin- 
Nernst signal from the bulk states through the Mott-like 
relation Therefore, in the next two subsections we 
will focus on analytical models to describe the scaling of 
the spin current and spin Hall conductance with the band 
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FIG. 6. Local spin current J z for normal metallic regime, 
and electrical bias, fi a — fi = — Qub — ft) with /i = 0.4io and 
T = 0. 



structure parameters and compare our intuitive analyti- 
cal models with the numerics. 

First however, to visualize the formation of the spin 
Hall effect at the sample boundary, it is instructive to 
plot the local spin current density in the numerical 4- 
band model. In order to do so, we first define a local 
spin current operator by 

jV) = -{[f,tfU(f-r')K, (12) 

i 

where u z is a Pauli matrix that acts on the spin space of 
the 4-band model. On the basis of the NEGF, it is then 
straightforward to evaluate the expectation value of the 
spin current operator at T = 0, which can be expressed 
as 



JZ ( r ) =J2 Tr [^'i^G^pG' 



(13) 



In Fig. [6j the local spin current is shown for the normal 
metallic bulk regime and an electrical bias applied from 
left to right. One clearly recognizes a spin current flowing 
along the edges. 

Note that our model shows local spin currents already 
at equilibrium. However at equilibrium, the spin current 
integrated over the cross section of a lead cancels and 
thus does not enter the spin Hall signal. For clarity, the 
local spin current that we show in Fig. [6] is only the 
non-equilibrium part. 

The rest of this section is organized as follows: parts A 
and B deliver two complementary approaches to explain 
the interface spin current that is transverse to the poten- 
tial gradient. In part C, the connection between these in- 
terface spin currents and the spin-Nernst coefficient gives 
us a qualitative understanding of the behaviour of spin- 
thermo effects for the bulk states. 



A. Effective 2-band models 

Here, we show that an anomalous spin-dependent ve- 
locity naturally appears within effective 2-band spin- 
diagonal electron or hole band models obtained by per- 
turbatively folding down the 4-band model (see Eq. IT])). 
We apply third order quasi-degenerate perturbation the- 
ory similar to our previous work [49] . 

The diagonal part of the Hamiltonian is Hq = 
diag(M, — M). For the perturbation part, we consider 

fct = ( A ~ + \ and = h'^*{~k), where V e/h = 

Ve/h T Bk 2 — Dk 2 and V e /h = Vq ± V m ; the subscripts 
e and h refer to electron and heavy hole bands, respec- 
tively. The B and D parameters will not enter to the spin 
current in third order perturbation theory. Note that, as 
compared to Eq. ([TJ , we allow for a finite in- plane poten- 
tial Vo in this analysis. Treating k± as operators acting 
on a perturbing potential, we obtain the spin-dependent 
effective 2-band Hamiltonians as follows (showing only 
the third order) 



= — 

eff,e gM 2 



2A 2 k+V h k_ -{A 2 k 2 ,V e }j , (14) 
2A 2 k_V e k+-{A 2 k 2 ,V h }) , (15) 



and ^off.o/h 



(^effc/h)*' ^ ne l° wes t order spin- 
dependent term of the effective electron/heavy hole 
model is thus given by 



fr-Pauli.c/h - ^ ( ft lff,e/h ~ h ttf,c/h 1 ° 



A 2 

'AM 2 



(V(±V -V m ) xk) z a z . (16) 



In the Heisenberg picture, we obtain a spin-dependent 
anomalous velocity 



an,e/h 



1 



v = Tr[r, hp au n e /h] 
in 



(17) 



A 2 ( 

Since \ an > e l h _|_ V(Vo =F V m ), we expect to see a spin 
current along the edge of the sample, similarly to the spin 
current carried by the edge states, but now the effect is 
induced by the bulk modes. 

Note that the assumptions for a valid perturbation the- 
ory are quite restrictive. The condition Ak <C 2\M\ re- 
stricts the energy range to about \E\ < 0.2to- Further, 
this approach works in the inverted regime only when one 
considers the bulk states and assumes a direct gap. The 
main drawback of this perturbative approach is, how- 
ever, that it assumes the variation of the potentials Vo, 
V m small compared to the gap 2\M\, which is not the 
case for the numerical confinement potential. Therefore, 
although we expect to find qualitative results by this 
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approach, it is important to compare it with the non- 
perturbative model including hard wall boundary condi- 
tions which will be done in the next subsection. 



B. Hard wall boundary spin current 

In this subsection, we present a complementary expla- 
nation of the spin current carried by the bulk states, valid 
also beyond the parameter regime Ak <C 2|M|, demon- 
strating that the reflection of an incident wave at a hard 
wall boundary leads to a spin current along the bound- 
ary. Due to a phase offset, this spin current persists even 
for a superposition of waves incident at different angles. 
We will show below that in the regime Ak -C 2M the spin 
current scales like A 2 /M . Interestingly, we observe that 
the explanation of the spin current given here seems to be 
close to what is seen in the numerical 4-band tight bind- 
ing model, because the numerically calculated spin Hall 
effect indeed scales like A 2 /M in the parameter regime 
Ak < 2M (with M > 0). 

In the model we consider now, the hard wall boundary 
condition for the envelope function is given by ip(y = 0) = 
0. While the direction of the outgoing beam is restricted 
by the energy and momentum conservation laws and is 
not spin-dependent, there is a spin-dependent phase shift 
between incident and reflected wave. Remarkably, even 
for an incident wave normal to the interface, a spin cur- 
rent moving along the interface is generated. In case this 
interface is bent, like it happens at the sample boundary 
connecting two perpendicular leads, it will transform into 
a spin Hall signal (like in Fig. ([6])). 

We start with the following ansatz for the spin f wave 
function 



As [eV nm/ft] 



4(y) = (v\tf) 



e ik « v u K 



(18) 



where the plane wave dependency on x has been sep- 
arated off. u± ky and u^x denote the spinors for fixed 
energy E and momentum k x . The condition ip^(0) = 
gives two equations for the coefficients r and c. The cor- 
responding coefficients for spin down can be found by re- 
placing k x — > ~k x . The operators of transverse velocity, 
V x(kx) = |fg and Vj(k x ) = -V2{-k x ) are indepen- 
dent of k y and complex-valued matrices. In the following, 
we will plot both spin up (in blue) and spin down currents 
(in red), evaluated by 

jlM = {^\5{y-y)V^), j^{y) = -jt fe ». (19) 
It is easy to see that the spin current 

f(y) = jI (y) - it (») = jI (y) + £ fc . (y) (20) 

is symmetric in the angle of incidence 9 = tan -1 k f L . 

Figure [7j shows the spin up and down currents for typ- 
ical parameters and the energy in the conduction band. 
The superposition of incoming and reflected propagating 




y [nm] 



FIG. 7. 



(Color online) The spin-dependent currents parallel 
to the interface at y — are plotted, with jl (y) in blue 
(solid) and (j/) in red (dashed). For perpendicular angle of 
incidence (k x = 0), the phase difference A(j} is always n. The 
energy is E = 0.3to in all cases. The thick/thin lines show 
the normal/inverted regime with M = O.lto and M — — O.lio, 
respectively. The other parameters are the same as for the 
lattice model. 



waves leads to an oscillating pattern. We are interested 
in the phase shift between spin up and down. The direct 
terms in jt (y) (i.e. two incoming or two outgoing propa- 
gating modes) are constant in y and current conservation 
dictates that the incoming and reflected currents are the 
same. Rotational invariance of the BHZ Hamiltonian and 
current conservation dictate that \r 2 \ = 1 independent of 
the spin. Because of time reversal symmetry, the current 
of the direct terms is independent of the spin and thus, 
the direct terms do not contribute to j s (y). 

The interference term between the incoming and out- 
going modes in jt (y) is given by 

2Re[(u kv \V x \u- ky )re~ 2ik y y ] 

= \(u ky \V x \u_ ky )r\2cos(2k y y - 4 - 4), (21) 



where <f>\ = a,rg((u ky \V x \u- ky )) and = arg(r). In Ref. 

[57] , 4>2~ 02 i s called the angle of giant spin rotation. At 
k x = 0, we have 



(u ky \V x \u- ky )\ K=Q 



A 2 kl + {M-Bk 2 ) 2 



(22) 



where k y is fixed by the energy. A first-order expansion 
in k, valid in the regime Ak <C 2\M\ yields 



js(y) oc 



A 2 k 



(23) 



In contrast to [58] , the spin current in our analytical anal- 
ysis is connected only with the propagating solutions as 
explained above. As one can see from Fig. [7J where the 
evanescent modes are included, the periodicity of j k (y) 
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FIG. 8. (Color online) The spin up phase shift of the re- 
flected current is <j>\ + 4>\- a) The phase (of the velocity 
matrix element) <j)\ as a function of the angle of incidence. 
It has the symmetry (p[(— 6) = it — 4>\(9). Different col- 
ors (blue, red, orange, ...) correspond to different choices of 
A = 0.05, 0.1, 0.4eV nm, where the limit A — > gives a 
step function. Parameter values are E = 0.3£o, M = O.lto 
and B, D have values as in Section II A. The picture does not 
change qualitatively, if we change the parameters of the un- 
derlying model, b) The phase (of the reflection coefficient) 4>X 
for the same parameters. For low values of A, it is propor- 
tional to A 2 , whereas for large values of A it saturates. 



and j\. (y) is only slightly affected which means that the 
evanescent contribution at least for the normal regime is 
minor and Eq. ( p3| still holds. 

Let us now analyze the phase relations between spin 
up and spin down currents more closely. The two phases 
<p[ and ipl (defined above) behave differently as a func- 
tion of the angle of incidence, as shown in Fig. [HI We 



find the symmetries <j)[{9) = n — 9) and <f> 2 {9) = 
-4>i(-0). For A -> 0, #(<?) = <p\{6) + becomes 
a step function, with <^(0) = ir/2. We are interested 
in A(f> — (fA — (f)+. For this, we again use a symme- 
try. If we flip the spin, r^(k x ) = r^(—k x ) implies that 
= 4>l(9) and V i {k x ) = ~V^(-k x ) implies that 



4(- 

i>{{-6) = TT + 4>\{e). Thus. 

A(j) 



+ = 2(0} + 01 ) (24) 



with A(j)(9) = 2ir — A<f>(—0). For not too small parame- 
ters A and small 9, the constant phase shift A<fi(0) = n 
is dominant. This phase shift ensures that the sign of 
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FIG. 9. (Color online) (a) Spin Nernst transmission function 
Tsn(^) for the normal regime (M — O.lto, blue dashed) and 
inverted regime (M = — O.lto, red solid). In the bulk gap, the 
edge states give rise to a quantized spin Hall conductance of 
e/2ir, except for the mini-gap. The black arrows are scaled 
according to the expected factor of ~™ h ~ 4.7 between con- 
duction and valence band spin Hall effect (see text), (b) Spin 
Nernst effect at T = 5 K, for normal (blue dashed) and in- 
verted regime (red solid). The latter is flipped horizontally, 
so we can compare signals of the same band character. 



the spin current is well-defined over a large range of 
9. Therefore, even the superposition of many incident 
modes at different angles (not included in this simple 
analysis) would lead to a well-defined sign of the spin 
current near the interface, while far from the interface, 
the spin current will be suppressed by the oscillations. 



C. Spin-Nernst signal for the bulk metallic regime 

In the preceding sections, we showed that the spin cur- 
rent can be understood by an anomalous velocity or a 
spin-dependent phase shift. The expressions we have ob- 
tained do not depend on the effective band mass (consid- 
ering the lowest order in Ak/M). We will now show that 
such a scaling of the spin current leads to a Nernst signal 
proportional to the effective band mass. 

Let us assume that the applied difference in the chem- 
ical potential A[i generates the spin and the charge re- 
sponses in the system. Then G s h — ^fG xx . For a given 
number of modes G xx is approximately constant. There- 
fore, using Eq. (23 1, one can see that G s h ~ m e / h A 2 /M, 
where m e /^ is the mass of the electron/ heavy hole band, 
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respectively. To the lowest order in Ak/M, the effec- 
tive 2-band and 4-band masses coincide and the 4-band 
effective masses m e / h are given by h 2 /2m e = —D — B 
and h 2 /2m,h = —D + B. Correspondingly, through 
the relation between the spin Hall conductance and the 
spin Nernst transmission signal Tsn ~ G s h ~ m e/h- 
The last dependence can be easily seen in the limit for 
D = 0, when the band structure of the BHZ model is 
particle-hole symmetric. Then, m e = —mh, which is 
consistent with Tsn ~ m e/h and the symmetry relation 
7sjv(m) — Tsn{~^) m that case. Figure [9] (a) shows 
numerical results for 7s at as a function of the chemical 
potential. In the 4-band model, the ratio of valence and 
conduction band effective masses is w 4.7. The 

black arrows are drawn for comparison of Tsn m con- 
duction and valence band and are scaled by the factor 
—mh/m e . Their position is chosen for energies corre- 
sponding to 4 propagating modes in the leads (counting 
spin), not counting edge states. For the normal regime 
(see dashed lines in Fig. [9] (a)) the scaling of the nu- 
merical Tsn is very close to what we predicted from the 
analytical approaches. In Fig. [91b), we show the corre- 
sponding spin Nernst signal. In the normal regime, we 
qualitatively find N s oc \nii\ (i = e,h) as expected from 
the Mott-like relation in combination with Fig. |9ja). 

In the inverted regime (solid red line) we must con- 
sider that near the bulk gap, the band character (E/H) 
has changed (compare the red/blue coloring in Fig. [2J; 
therefore, the band for \i > gets a heavy hole character. 
Further, as long as the edge states do not yet merge to the 
bulk, they are responsible for an offset of Tsn = —2. The 
black arrows again indicate the factor that we ex- 
pect for the comparison of conduction and valence band 
signals at the same number of contributing modes, how- 
ever now we are measuring the signal from the level of 
the edge states. Analyzing numerically the scattering 
matrix we find that in the valence band the contribution 
to Ts n 01 bulk and edge states are additive, while this 
is not the case for the conduction band. Taking into ac- 
count this fact, it is suprising that the simple analytical 
analysis applicable to the normal regime still describes 
qualitatively the numerics. We believe that this might be 
the case, because the first bulk state resembles the edge 
state character, and our argument about the symmetry 
of Tsn(ij-) — Tsn{—^) for the particle-hole symmetric 
Hamiltonian still holds. 



V. CONCLUSION 

We have analyzed the thermoelectric transport in four- 
terminal setups of HgTe/CdTe quantum wells with a par- 
ticular emphasis on spin-dependent effects due to spin- 
orbit coupling. Thereby, we have used a combination 
of analytical and numerical methods to analyze spin- 
dependent transport phenomena. The Seebeck and the 
spin Nernst signal show a peculiar dependence on the pa- 
rameters of the Bernevig-Hughes-Zhang model which can 



be qualitatively understood as originating from a spin 
Hall effect that arises at in-plane potential or confine- 
ment boundaries of the system. We have demonstrated 
that the spin Nernst effect is a strong experimental tool to 
get a better understanding of the mini-gaps that arise due 
to the spatial overlap of edge states on opposite sample 
boundaries. Most interestingly, we have derived a Mott- 
like relation between the spin Nernst coefficient and a 
smoothed spin Nernst transmission function that is valid 
to all orders in the temperature difference between the 
warm and the cold reservoir. Our findings might help to 
optimize future experiments on thermoelectric transport 
properties of two-dimensional topological insulators. 

Financial support by the German Science Foundation 
(DFG, SPP 1285), the Helmholtz Foundation (VITI), 
and the European Science Foundation (ESF) is gratefully 
acknowledged. 



Appendix A: Tight binding Hamiltonian 

Using the representation of the plane wave annihila- 
tion operator in the basis of lattice sites (J, I), Ck = 



^3- J2j 1 e*°( *^ + * !t ' * 'Cjj, we obtain the following substitu- 
tion rules for the continuum model momentum operators, 



1 



d k — sin(afe x )clck 
bz a 



\_ \ " t i 



2ia 



j.i 



k l / d 2 k - (2 - 2 cos(afcz)) c^Ck 



(Al) 



BZ " 

; ~a? c 5+i,l c W ~ 2c h c J,l + c ]-i,l c 3,U (A2) 
U 

and analogous rules for k y and k 2 . 

We define matrices on the band space T 1 = r x a z , T 2 = 
-T y a , T 3 = r -^a v T 4 = -^±^a x and T 5 = T z a 
where Tj are the Pauli matrices acting on the E/H space 
and <Ji acts on the spin space (to,cto are unit matrices). 
Then, the lattice Hamiltonian corresponding to H of Eq. 
([I]) reads 



tb 



^ r [AT 1 +R Q T 3 BT 5 + D1\ t 

2^ . — + -—2 — JW*« 



2ia a 2 

AT 1 + i? r 3 BT 5 + Dl 

1 . 

2ia a 2 

AT 2 + R T 4 BT 5 + Dl 



2ia 



AT 2 + R T 4 BT 5 +D1\ t 

C j,l-l C J,l 



2ia 



r(BT 5 + Dl) + MT 5 ) c^Cjj 



(A3) 
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where the summation over the grid points (j, I) is re- 
stricted by the geometry of the sample and we only in- 
cluded the most important linear Rashba term, propor- 
tional to i?o m Eq. 



Appendix B: Mott-like relation 

In this Appendix, we show how the Mott-like relation 
of Eq. |9) can be generalized to finite temperatures. 
For this we consider the Fourier transform of transmis- 
sion functions AT c ,„(t) = i / dEe~ iEr AT cq . The spin 
Nernst effect is defined as 



N s = 



2AT 8ttAT 



at/ 



d£^AT c ,„(£)(/ c 



ft), (Bl) 



where the potential /i is assumed to be the same for all 
leads, while the temperatures may differ. The integral 
has the form of a convolution. The Fourier representation 
is 



N s = 



1 

4tt 



dr—e 



AT c ,,(r) 



x 



^ AT sinh x 
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x—irrkB^c 

X — TTTkB^o 



(B2) 



Further, we define a symmetric "smoothing" function that 
depends on temperatures of the leads c and q as 



F q {r) 



7r 2 fc|T c AT sinh a; 



X — TTTkB^o 



(B3) 



If we put T c - T, = AT, we find with AT ^ 



F q (r) « 3 



x coth x — 1 



■ sinh ; 



(B4) 



=7rrfc B T c 



which has a width of At sa , 4 „ and F q (0) — > 1. Now, 

7T Kb -I- c V ' 

we define a temperature-smoothed spin Nernst transmis- 
sion function as 



fs N (T)=J2F q (r)AT C}q (r), 



(B5) 



which implies that Tsn(E) = f dre iEt Tsn{t) is real. 
Finally, we obtain the relation 



vrfc|T c df SN {E) 
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dE 



(B6) 



which is exact to all orders in T c and AT. The meaning 
of the latter equation is the following one: First taking 
the derivative OeTsn and then smoothing by temper- 
ature is the same as first smoothing with a modified 
smoothing kernel and then taking the derivative. 
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